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Abstract 

We present results of a study of the charmless vector-vector decays —>■ K*^p^ and B~^ — > 
K*^p+. The results are based on a 140 fb ^ data sample collected by the Belle detector at the 
KEKB asymmetric e'^e" collider. We obtain the branching fraction B{B~^ — ^ K*^p'^) = (6.6 it 
2.2(stat.) ± 0.8(syst.)) x 10~^, and set upper limits on the branching fractions B{B'^ K*^ p^) < 
2.6 X lO"*' and B{B^ fo{980)K*^) < 5.2 x 10"^. We also perform a helicity analysis of the p and 
K* vector mesons in the decay B~^ K*^p'^, and obtain the longitudinal polarization fraction 
Ro{B+ K*^p+) = 0.50 ± 0.19(stat.)l[j:[Jf (syst.). 

PACS numbers: 13.25.Hw, 14.40.Nd. 
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In addition to rate asymmetries, B — > VV decays provide opportunities to search for 
direct CP and/or T violation fi\ and new physics i.2] through angular correlations between 
the vector meson decay final states. These decays produce final states where three helicity 
states are possible, whose amplitudes are called Aq, and A_ in the helicity basis. The 
Standard Model (SM) with factorization predicts > 0,3, where i?o = |^olV(l^oP + 
|yl+|2+ {Rt = (|yl+p + |y4_p)/(|AoP + |A+|2+|A_p)) is the longitudinal (transverse) 
polarization fraction. In the tree-dominated — p^p°, this prediction is confirmed 0,1^. 
In contrast, for the B (pK* decay, which is a pure b s penguin transition. Belle |7| 
and Babar [6] find Rq ~ Rt, which is in disagreement with SM predictions. It is thus 
important to obtain polarization measurements in other VV modes, such as -B ^ K*p and, 
in particular, in the pure penguin b — >■ sdd decay, B~^ — *■ K*^p'^ . 

In this paper, we present the results of a study of B^ — »• K*^p^ and B^ K*^p^ decays 
with a 140 fb~^ data sample containing 152 x 10^ B meson pairs collected with the Belle 
detector at the KEKB asymmetric-energy e~^e~ collider [s"] operating at the T(45') resonance 
{^/s = 10.58 GeV). The production rates for B^B~ and B^B^ pairs are assumed to be equal. 

The Belle detector is a large solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in detail elsewhere j^. 

We select B K*p candidate events by combining four charged tracks (three pions plus 
one kaon) or three charged tracks (two pions plus one kaon) and one neutral pion. Each 
charged track is required to have a transverse momentum px > 0.1 GeV/c and to have 
an origin within 0.2 cm in the radial direction and 5 cm along the beam direction of the 
interaction point (IP). 

Particle identification likelihoods for the pion and kaon particle hypotheses are calculated 
by combining information from the TOF and ACC systems with dE/dx measurements in the 
CDC. To identify kaons, we require the kaon likelihood ratio, Lk/{Lk + L^^), to be greater 
than 0.6. To identify pions, we require Lk/{Lk + -^tt) to be less than 0.4. In addition, 
charged tracks are rejected if they are consistent with the electron hypothesis. 

Candidate 7r° mesons are reconstructed from pairs of photons that have an invariant 
mass in the range 0.1178 — 0.1502 GeV/c^, corresponding to a window of iScr around the 
nominal 7r° mass, where the photons are assumed to originate from the IP. The energy of 
each photon in the laboratory frame is required to be greater than 50 MeV for the ECL 
barrel region (32° < 6 < 129°) and 100 MeV for the ECL endcap regions (17° < < 32° 
or 129° < 9 < 150°), where 9 denotes the polar angle of the photon with respect to the 
beam line. The 7r° candidates are kinematically constrained to the nominal vr*^ mass. In 
order to reduce the combinatorial background, we only accept vr^ candidates with momenta 
Pt,o > 0.40 GeV/c in the e~^e~ center-of-mass system (CMS). 

Candidate p mesons are reconstructed via their p° — >■ tt+tt" and p+ tt+tt^ decays. For 
both the charged and neutral modes, we require 0.62 GeV/c^ < M^nn) < 0.92 GeV/c^. We 
select K*^ — >■ K~^tt~ decay candidates with invariant masses in the range 0.84 GeV/c^ < 
M{K+TT-) < 0.94 GeV/c2 for /s:*V° and 0.83 GeV/c^ < M{K+n-) < 0.97 GeV/c^ for 
K*'^p+. 

To isolate the signal, we form the beam-constrained mass Mbc = y E^^^^ — Pb-, and the 
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energy difference AE = Eb — -Ebeam, where -Ebeam is the CMS beam energy, and pb and Eb 
are the CMS momentum and energy, respectively, of the B candidate. We accept events 
in the region defined by Mbc > 5.2 GeV/c^ and -0.4 GeV < AE < 0.4 GeV. Within this 
accepted range we further define mode-dependent signal regions. For K*^p^, the AE 

signal region is -0.04 GeV < AE < 0.04 GeV. For B+ K*^p+, the AE distribution has 
a tail on the lower side caused by incomplete longitudinal containment of electromagnetic 
showers in the CsI(Tl) crystals, so the AE signal region is broadened to —0.10 GeV < AE < 
0.06 GeV. For both decays, the Mbc signal region is 5.27 GeV/c^ < Mbc < 5.29 GeV/c^. 
These requirements correspond to approximately ±3cr for both quantities. 

The continuum process e~^e~ qq {q = u, d, s, c) is the main source of background and 
must be strongly suppressed. One method of discriminating the signal from background 
is based on the event topology, which tends to be isotropic for BB events and jet-like 
for qq events. Another discriminating characteristic is 6b, the CMS polar angle of the B 
flight direction. B mesons are produced with a 1 — cos^ 6b distribution while continuum 
background events tend to be uniform in cos6b- For B^ —>■ K*^p^, we require | cos^thrl < 0.8, 
where ^thr is the angle between the thrust axis of the candidate tracks and that of the 
remaining tracks in the event. This distribution is fiat for signal events and peaked at 
cos6'thr = ±1 for continuum background. 

We use Monte Carlo (MC) simulated signal and continuum events to form a Fisher dis- 
criminant based on modified Fox- Wolfram moments that are verified to be uncorrelated 
with Mbc, AE and variables considered later in the analysis. Probability density functions 
(PDFs) derived from the Fisher discriminant and the cos6b distributions are multiplied 
to form likelihood functions for signal {Cg) and continuum (Cqq); these are combined into 
a likelihood ratio TZg = Cs/{Cs + Cqq). Additional discrimination is provided by the h- 
fiavor tagging parameter r, which ranges from to 1 and is a measure of the likelihood 
that the h flavor of the accompanying B meson is correctly assigned by the Belle flavor- 



tagging algorithm Events with high values of r are well-tagged and are less likely 

to originate from continuum production. We define a multi-dimensional likelihood ratio 
M = £f ^^^/(/:f + where C^/^^^ denotes the likelihood determined by the 

r-TZs distribution for signal and Cq^^^^ is that for the continuum background. We determine 

A4 cut by optimizing the figure of merit, S/ ^{S + B), where B is the number of background 
events and 5" is the number of signal events. We require Ai > 0.85 for K*^p^ and M. > 0.95 
for K*°p+. 

To investigate backgrounds from other B decays, we use a sample of BB MC events 
corresponding to an integrated luminosity of 412 fb~^. We find background from B~^ 
D^{K*^TT^)7r'^ in the p sideband region. We apply the requirement |M(i^'7r"'"7r'') — Mo"\ > 
0.050 GeV/c^ to veto those events. This requirement does not remove any B^ K*^ 
events. The B^ K*^p^ mode includes vetoes on B^ (K^ 71^71^)71^ , B^ — >■ 

L'*(2010)-(L'°7r-)7r+ and D*(2460)~(L'°7r")7r+. The cuts \M{K7nr) - Mdo\ > 0.050 

GeV/c2, \M{K7ni) -Md" {2ow) \ > 0.050 GeV/c^ and \M (Ktiti) - Md* (^2460) \ > 0.060 GeV/c^ 
eliminate these backgrounds. 

After all the event selection requirements are applied, there is no significant peaking 
background in either AE or Mbc foi' —>■ K*^p^ . For B^ K*^p^ we find some broadly 
peaking 3-body and 5-body rare decays above and below the signal in the AE distribution. 
These shapes are included in the M\^c-AE fit for this mode. A further cut, AE > —0.2 GeV, 
is also applied to eliminate most of the remaining 5-body background. 

In the B^ K*^p^ A/bc-Ai? signal region, we find that the fraction of multiple candi- 
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dates is 3.5% for the \Aq\'^ helicity state and 1.6% for the |A±P state. We allow for multiple 
candidates in this mode. For K*^p^ the wide mass windows result in many multiple 

candidates; we choose the entry that has the minimum from the vertex constrained fit. 

We extract the signal yield by applying an extended unbinned maximum-likelihood fit 
to the two-dimensional Mbc-A£' distribution. The fit includes components for signal plus 
backgrounds from continuum events and 6 — c decays. The signal PDF is represented by a 
Gaussian function for Mbc and either a double Gaussian (for the K*^p^ mode) or a "Crystal 
Ball" line shape function (for the K*^ mode) for ^E. The shape parameters are 
determined from fits to MC. The signal PDF is adjusted to account for small differences 
observed between data and MC. The — > K*^p^ signal PDF is calibrated with B^ — >■ 
D^[K~n~^7[~)7!'~^] the B^ — > K*^p'^ PDF is calibrated with high-statistics modes containing 
7r° mesons, i.e., B^ — > D^{K~^h~ti^}%~^ . The continuum PDF is described by a product 
of a threshold (ARGUS) function [1^ for Mbc and a first-order polynomial for AE, with 
shape parameters allowed to vary. The PDF for b ^ c decay is modeled by a smoothed 
two-dimensional histogram obtained from a large MC sample. In addition, the B^ K*^p^ 
fit contains PDFs for the dominant rare 3-body {B^ — > K*^tt) and 5-body {B^ — ^ aiK*^) 
decays, modeled by smoothed two-dimensional histograms. 

Figure^lshows the fit results. We find 14.1 ±4.4 B^ K'^Ti^n^n^ events and 56.5 ±11.6 
B~^ — >■ K^n'TT^TT^ events. The corresponding statistical significance of the signal, defined 

as ^—2 ln(£o/'^max), where £max is the likelihood value at the best-fit signal yield and Cq 
is the value with the signal yield set to zero, is 5. Oct or 6. Sex. 

The Mbc ~ fits do not distinguish signal from non-resonant decays, such as B ^ pKir 
OT B ^ K*TT7T, which contain the same final state particles. Therefore, we extract the K*p 
signal yield by fitting the M{7nr) and M{Kti) invariant masses for events in the Mbc and 
AE signal region. 

For B^ — >■ K*^p^ signal region events, we perform a two-dimensional maximum likelihood 
fit in the range 0.5 GeV/c^ < M{K+n-) < 1.2 GeV/c^ and 0.5 GeV/c^ < M(7r+7r-) < 
1.2 GeV/c^. The distribution for these events is shown in Fig. |21 From this distribution it 
is clear there are significant contributions from non-resonant backgrounds. The signal PDF 
is modeled by the product of Breit-Wigner functions for both the p° and K*^ resonances 
with means and widths determined from the MC simulation. The continuum and b — >■ 
c fractions are each modeled by a smoothed two-dimensional histogram determined from 
MC. Three main non-resonant decays are considered in this fit: B^ K*^7t~^7t~, B^ — »• 
Ki{imy{K*^n+)7i- and B° p°K~n+. The K*°7i+n- and p^K~'n+ components are 
represented as a Breit-Wigner function for the resonance, multiplied by a second order 
polynomial fit to the non-resonant phase space, determined from MC. The i^'i(1400)"'"7r~ 
component is modeled by a smoothed histogram. From the data mass distribution it is 
apparent that an additional component to model the B^ —>■ /o(980)i^'*° decay is needed. 
This is represented by the product of two Breit-Wigner functions for the K*^ and /o(980) 
resonance. The K*^ mean and width are determined from MC. The /o(980) mean and width 
are allowed to fioat within the range of values reported in the PDG jl^. The continuum 
and 6 — s> c backgrounds are fixed at the levels determined from the Mbc - fit; the other 
components are allowed to fioat. The projections of the fit are shown in Fig. El Table|l]lists 
the results of the fit. The statistical significance of the fo{980)K*^ signal is 2.8a. 

For B^ K*^p^, we fit M(7r7r) in the K* signal region and M{Kti) in the p signal 
region simultaneously. Figure El (bottom) shows the fit result. The solid curve shows the 
result of a binned maximum-likelihood fit with four components: signal, background from 
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FIG. 1: Projections of Mbc for events in the AE signal region (left), and projection of AE in the 
Mbc signal region (right). The upper plots are for — > K*^p^, the lower plots for K*^p^ . 

The solid curves show the results of the fits. The dotted lines represent the continuum background. 
The sum of 6 — > c and continuum background component is shown as dashed lines. 



continuum and BB^ and non- resonant pKi: and K*^7nT. The signal p and K*^ components 
are represented by Breit-Wigner functions with masses and widths determined from MC 
simulation. The background from continuum and BB is described by a threshold function 
plus a Breit-Wigner function for M(7r+7r°), a threshold function plus a Gaussian and a 
Breit-Wigner function for M(i^'+7r~), where a resonant component is included to account 
for resonance production in the continuum. The shape parameters are determined from 
sideband data. The non-resonant pKn and K*7Tn components are represented by threshold 
functions with parameters determined from a MC simulation with the final state particles 
distributed uniformly over phase space. In the fit, the yield for each component in the 
p or K*^ signal region is required to be the same. All normalizations are allowed to float, 
except for the background from continuum and BB, which is fixed at the Mbc-Ai? fit results. 
Table ITTl lists the yields in the p~^ and mass window. The statistical significance of the 
B+ K*^p+ signal is 3.2a. 

We use the p'^ 7r+7r° and K*^ K~^7i~ helicity-angle (6'hei) distributions to determine 
the relative strengths of \Aq\'^ and |A±p. Here 6'hei is the angle between an axis anti-parallel 
to the B flight direction and the vr"*" {K~^) flight direction in the p"'' (K*) rest frame. For the 
longitudinal polarization case, the distribution is proportional to cos^ Ohei{p) cos^ Ohei(K*) and. 
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M(7tV) (GeV/c^) M(KV) (GeV/c^) 

FIG. 2: Distributions of M(7r+7r~)(left) and M{K~^tt~) (right) in the K*^p'^ signal region (open 
histogram) and Mbc sideband (hatched histogram). Peaks for and K*^ can clearly be seen and 
a sharp peak for the /o(980) is apparent. Significant contributions of non-resonant background are 
also seen. 



TABLE I: Results of fit to invariant mass for B K p . The yields are the number of events 
over the whole M{'ktt) - M(Ktt) region. 





/o(980)/s:*° pKtt 


K*°Tnr KiTT |yloP effic.(%) 


|vl±|2 effic.(%) 


0±5.2 


10.2t|:3 30.5^1^ 


22At\li 9Atli 1.9 


3.4 



for the transverse polarization case, is proportional to sin^ 6'hei{p) sin^ 6'hei(i^*), where ^hei(p) 
(6'hei(_ft'*)) is the helicity angle for p {K*). Figure HI shows the distributions of the cosine of 
the helicity angle for p and for K* for events in the Mbc and AE signal region. We perform 
a two-dimensional unbinned maximum likelihood fit to the p and K* helicity distributions. 
The fit includes components for signal, backgrounds from continuum and BB and non- 
resonant pKn. The K*'ktt component obtained from the invariant mass fit is small and is 
not included in the helicity fit. PDFs for signal \Aq\'^, |^±P helicity states are determined 
from the MC simulation. The PDF for background from continuum and BB is obtained 

TABLE IL Results of fit to invariant mass for K*^p^, together with the MC-determined 

efficiencies for |Aop and |A±P states. The yields are the number of events in the K* and p mass 
window. 





pKir 


K*-K-K 


|Aop effic.(%) 


|A±|2 effic.(%) 


26.6 ±8.7 


12.8 ±3.3 


0.6 ±3.7 


2.0 


3.3 
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M(7i^7i") (GeV/c^) M(KV) (GeV/c^) 




M(7r^ tt" ) (GeV/c ) M(K+ n ) (GeV/c) 



FIG. 3: The upper plots are for — > K*^p^, M(7r^7r^) (left) in the K*^ signal region and 
M{K~^'iT^) in the p signal region (right). Solid curves show fit results. Dotted lines show continuum 
and non-resonant components, red dashed line shows the /o(980)i^*'' component. The lower plots 
are for — > K*^p~^ , M{'k~^t:^) (left) in the K*^ signal region and M{K~^'k~) in the p signal 
region (right). The solid curves show the results of the fit. The signal (continuum, pKn, K*^7nT ) 
component is shown as a red dashed (dot-dashed, dotted, dashed) line. 

from sideband data. The pKn PDF is determined by fitting the hehcity distribution for 
events in the region 1.0 GeV/c^ < M{K7i) < 1.5 GeV/c~^; it is consistent with a cos^ 6'hei- 
hke cos6{p) hehcity and a flat cos^(i^7r) distribution. The fit result, shown as the solid 
histogram, gives Rq = 0.50 ± 0.19; 3.2cr away from 100% longitudinal polarization. 
We obtain the longitudinal polarization fraction, 

Ro{B+ K*°p-^) = 0.50 ± 0.19(stat.)^[J:|]^(syst.), 

where the systematic error is associated with the fitting procedure determined by shifting 
each parameter by ±lcr, changing the PDF, and taking the quadratic sum of the resulting 
changes in Rq as the systematic error. 

We use the mass fit result and MC-determined efficiencies weighted by the measured 
polarization components to calculate the branching fraction for — > K*^ for the cal- 
culation of the upper limit of — > K*^p^, 100% longitudinal polarization is assumed. 
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FIG. 4: Projections to cos0hei(p) (left) and cos^hei(K*) (right). The solid histograms show the 
results of the two-dimensional fit. The red dashed (blue dot-dashed) histograms are the Aq {A±) 
component of the fit; The dotted histograms are backgrounds from continuum and BB. The black 
dashed histograms are for non-resonant pKn. The low event yield near cos 0hei(p) = 1 is due to the 
p^o > 0.4GeV/c requirement. 



We consider systematic errors in the branching fraction of the decay B K*p that 
are caused by uncertainties in the efficiencies of track finding, particle identification, 7r° 
reconstruction, continuum suppression and fitting. We assign a 1.1%/track error for the 
uncertainty in the tracking efficiency. This uncertainty is obtained from a study of partially 
reconstructed D* decays. We also assign a 0.6%/track error for the particle identification 
efficiency that is based on a study of kinematically selected D*~^ D^tt'^, K^tt^ 
decay. A 4.0% systematic error for the uncertainty in the 7r° detection efficiency is determined 
from data-MC comparisons of ^ Tr^vr'^Tr*' with 7] tt+tt^tt'^ and 7] —>■ 77. A 3.4% {K*^p^)/ 
3.8% {K*^p'^) systematic error for continuum suppression is estimated from studying — >■ 
D-7r+, D- K-ir+n- and B+ Wti^, W K+n-n'^. A 3.8% systematic error 
associated with the K*^p^ fit is obtained by shifting the parameters by ilcr and changing 
the PDF. A fit systematic error of 16.1% for is:*V° and 7.8% for /o(980)/s:*° is found by 
varying the fit parameters by ±la. A 0.5% error for the uncertainty in the number of BB 
events in the data sample is also included. For B^ —>■ K*^p^, we also include a 11.9% error 
due to the uncertainty in the fraction of longitudinal polarization. The quadratic sum of all 
of these errors is take as the total systematic error. We obtain the branching fraction 

B{B+ K*^p+) = (6.6 ± 2.2(stat.) ± 0.8(syst.)) x 10"^ 

and set 90% confidence level (CL.) upper limits on B^ K*^p^, 

B{B^ K*^p°) < 2.6 X 10-^ 

and 50 /o(980)ir*o, 

B{B^ fo{980)K*^) < 5.2 x 10"^ 

In summary, we measure the branching fraction for B^ —>■ K*^p^, set a 90% CL. upper 
limit on B^ K*'^p^ and B^ fo{980)K*^. A heficity analysis is performed for 5+ — »• 
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K*^p'^ . We find a substantial transversely polarized fraction. The results are consistent 
with some recent calculations 0, 0| . All results reported here are preliminary. 
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